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Abstract 

In this work we calculate the direct solution of the equations for the retarded 
electron-boson interaction in the case of (f-wave symmetry for the pair wave 
function and in the real axis formulation. We use a spectral function contain- 
ing an isotropic part and an anisotropic one: a 2 (Q, <ft, (j>')F(Q) = cQF(fl) + 
a d F(Q) • cos (20) cos (2<j)') and make the simple assumption: a d F(Q,) = 
gdalF(Q) where gd is a constant. For appropriate values of the isotropic 
electron-boson coupling constant \ s and the anisotropic one A^, solutions are 
obtained with only d-wave symmetry for the order parameter and only s-wave 
one for the renormalization function. We have employed the real axis formu- 
lation in order to compare the theoretical curves to the tunneling density 
of states of the optimally-doped high-T c superconductor F^SraCaC^Og+a;. 
The results of our numerical simulations are able to fit very well the value 
of the gap, the critical temperature and the shape of the density of states in 
all the energy range, as recently determined in our break-junction tunneling 
experiments. At T > T c the theoretical conductance still shows a broaden 

peak that disappears at T* ~ 140 K. 

PACS numbers: 74.50.+ r; 74.20.-z; 74.72.Hs 
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At the present moment, there are many experimental results that strongly suggest the 
presence of a dominant d-wave symmetry of the pair wave function in high-T c cuprates, at 
least the layered ones Recent tunneling measurements on the layered high-T c super- 
conductor Bi 2 Sr 2 CaCu208+:); (Bi 2212) yielded high or very high values of the low- 
temperature gap, gapless features and a small but finite value of the quasiparticle density of 
states (DOS) at zero bias. These facts could be explained, at least partially, by assuming a 
(i-wave pair state which has nodes within the superconducting gap. On the other hand, in 
the last years, a few attempts have been made to fit some of the high-T c superconducting 
properties by using the Migdal-Eliashberg theory P-|llfl, but most of the works concerned 
only the solution in s-wave pair symmetry. In particular, we recently used the supercon- 
ducting DOS of optimally-doped Bi 2212 single crystals with T c ~ 93 K, determined in 
break-j unction tunnel experiments, in order to reproducibly determine the electron-boson 
spectral function of Bi 2212 by inversion of the s-wave Eliashberg equations 0JH]. 

In the present paper, we extend the real-axis direct solution of the Eliashberg equations 
for strong electron-boson coupling to the d-wave case and discuss the main features of the 
solution as function of temperature and boson coupling strength. In the framework of 
this model, we are able to obtain a very good agreement between our experimental break- 
junction tunneling data and the theoretical results with regard to both the DOS shape in 
all the energy range and the critical temperature. 

Our numerical analysis starts from the well-known generalized Eliashberg equations for 
the order parameter A(u, k) and the renormalization function Z{u,\z). They have kernels 
that depend on the retarded interaction a 2 (Q, k, k.')F(Q), the coulomb interaction /i* (k, k'), 



and the effective band fl2|--|15|| . For simplicity we suppose that k and k' lie in the Cu02 
plane (ab plane) and we neglect the relatively small band dispersion and the gap in c direc- 
tion. Therefore we use a single band approximation where the Fermi line is nearly a circle 
while <p and (f)' are the azimuthal angles of k and k' in the ab plane, respectively. We thus 
expand a 2 (Q, <p, <p')F(Q) and //(</>, </>') in terms of basis functions ipi (0) where the first few 
functions of lowest order are: 



fa (0) = 1; fa (0) = cos (20); ^ 2 (0) = v/2 sin (20); 

In the present paper we study the simplest d-wave electron-boson model interaction 
expressed by the following separable expressions 



a 2 (ft, 0, <f>')F{n) = a 2 00 F(n) + a 2 n F{tt) ■ fa (0) Vi 



(1) 



At* (0,00 = / i* + yU * 1 • ^ (0)Vi (0') 



(2) 



and we search for a s + id solution given by 



A(w, 0) = A t (u) + A d {uj) ■ fa (0) 
Z(u, 0) = Z s (u) + Z d (u) ■ fa (0). 



With these assumptions the Eliashberg equations in presence of impurities [IE| become: 



AJuj)ZJuj) 



+oo 



P s (uj') ■ [K 00+ (uj, lo') - /Xqo • © (u c - u')] duj' + 



mr-P s (cu)/ C 2 + P s H + iV s (o;; 



r+oo 

A d (u)Z s (u) = / P d (u/) • [K 11+ {u,u') - /i* n • 6 (cu c - w')] da;' 



(3) 



(4) 



/■+oo 

[1 - Z s {u)\ -u= N 8 (u') ■ K m _{u,u')dJ + 
Jo 

iirT ■ N s (uj)/ \C 2 + P s (lo) + N s (uj) 



(5) 



where T = n I / 7V(0)7r 2 , C = cotg(8 ), rij is the impurity concentration, N(0) is the value 
of the normal DOS at the Fermi energy and Sq is the scattering phase shift. We have not 
written the equation for Z d (u) that is homogeneous. In the weak-coupling case its only 
solution is Z d {uS) = 0, while in the strong-coupling one, in principle, there is a possibility 
that a nonzero solution exists above some coupling-strength threshold. In the present work, 
we assume that the stable solution corresponds to Z d (u) = for all the couplings and we 
do not consider the rather exotic possibility of Z d {uS) ^ JIBUTI] . In eqs. (3)-(5), u c is a 
cut-off energy and we have 



r+oo 

K ii± (u,u')= / dnal(n)F(n)- 

Jo 

fi + a;' + a; + i • 5 VL + u' — uj — i-5 Q — uj' + uj + i ■ 5 Q — uj' — uj — i ■ 5 

where / (a/) and n (fl) are the Fermi and Bose functions, respectively, and furthermore we 
know that 

P s (u) = (l/2vr) ^ Real (a s {uj) / ^ uj 2 - A 2 (uj) - [A d (uj) ■ fa (0)] 2 ) d0 

PiH = (1/2tt) jf** Pea/ (A d (u;) • ^ (0) (0) / \j ^ - A 2 s {uj) - [A d (uj) ■ fa (0)] 2 ) d0. 

The quasiparticle DOS that is compared to the experimental data is given by 

N s (oj) = (1/2tt) Real L/sfu 2 - A*(w) - [A d (u) ■ fa (0)f ) d<j>. (6) 

In the case we want to obtain a s + d solution, we need only to replace the denominator of 
P d {u), P s (uj) and N s (uj) with \J uj 2 - [AJu) + A d (uj) ■ fa (0)] 2 . 

In our numerical analysis we have put, for simplicity, a 2 1 F(Q)=g d a 2 l0 F(Q) where 
g d is a constant |T7|]. As a consequence, the electron-boson coupling constants for 
the s-wave channel and the d-wave one are A s = 2 dVta^F{Q) /Vt and = 
(1/tt) jf* d<j>fa (0) J +o ° dna 2 n F{n)/n = g d ■ X s , respectively. 

We solved the generalized, real-axis, Eliashberg equations (3)- (5) in direct way by using 
an iterative procedure that continues until all the real and imaginary values of the functions 
A s (u), A d (u) and Z s {u) at a new iteration show differences less than 1 • 10~ 3 with respect 
to the values at the previous iteration. Usually the convergence occurs after a number of 
iterations between 10 and 15. By using this approach we found that A(u, 0) has either a 
pure 5-wave or a pure d-wave symmetry depending on the values of the coupling constants 
X s and X d . The values of the parameter g d = X d /X s determine the stability of the two types 
of solution: if 1.5 < g d < 5, there are no problems for an automatic convergence to the 
pure d-wave solution. When g d < 1.5, the choice of the starting values of A s (u) and A d (uj) 
determines the symmetry of the final solution to which the iteration procedure converges. 
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If gd 3> 1 (i.e. roughly > 7) the iteration procedure converges only in the complex-axis 
formulation. 

We can now check if exists a couple of X s and A^ values that, together with a proper 
electron-boson spectral function a 2 (Q, (ft, <p')F(Q), can reproduce our tunneling density of 
states N ex (u) of Bi 2212 in the framework of a pure c?-wave strong-coupling solution. For do- 
ing this we compare the N ex (uj) to the quantity N(u) = N s (E+lu) \f(E)-f(E + u)\dE 
where N s (u) is calculated from eq.(6) by using al F(Q) = (A s /A) a 2 E(Q)Bi22i2- As a first- 
order approach to the d-w&ve modeling of the tunneling curves, we used the electron-boson 
spectral function previously determined by the inverse solution of the s-wave Eliashberg 
equations applied to the same Bi 2212 data as a 2 E \Q) smxi ( see the inset of Fig. 1) ||f|. 
Of course, A is the corresponding coupling constant. 

In Figure 1 our tunneling experimental data (open circles) and the best fit curve at 4 K 
(solid line) that is obtained for A s = 2, = 0, gd = 1.15 and yields almost the experimental 
T c , T^ alc = 99 K, are shown. We want to emphasize the impressive agreement of the 
theoretical curve to the experimental one in all the energy range, including the region at V < 
A M (peak of the conductance) where states inside the gap in all the experimental tunneling 
curves are present. Taking into account a small amount of impurities corresponding to 
T = 0.3 meV and C ~ 0, also the zero-bias conductance is reproduced and T c caZc = 96 K. 
In the case of g& = 1.12 and no impurities, T c mk = 93 K, that is exactly the experimental 
value, but the fit is a little less good. It is interesting to note that there are different 
pairs of A s and X d values with 2 < A s < 2.5 and 2 < X d < 2.25, all giving T c caZc = 93 
K and producing a good fit of the experimental data, the best remaining that one shown 
in Fig.l. When the <i-wave component of the Coulomb pseudopotential is different from 
zero (the s-wave component does not affect the results in pure (i-wave symmetry) the above 
mentioned best-fit parameters give the experimental T c value for = 0.0125. The short 
dash line of Fig. 1 shows the normalized dl/dV curve obtained at 4 K in this case. It is 
practically indistinguishable from the curve at /i^ = 0. In general, the main consequence 
of a fi*i 7^ is the reduction of the values of T c caZc and Am- In Fig. 1 we also show the 
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theoretical normalized conductances, obtained by using the parameters of the best-fit curve, 
at different temperatures between 60 K (long dash line) and 140 K (tiny short dash one). 
We can see that the normalized conductance is practically flat only at T* ~ 140 K > T c : 
this is a strong-coupling effect that occurs at A > 2 and has been briefly discussed in s-wave 
symmetry JT^-^2[ . Here we demonstrate that it also exists in c?-wave symmetry, is enhanced 
by the particular characteristics of Bi 2212 and could be related to the pseudogap observed 
in underdoped and overdoped Bi 2212 ||. 

In this <i-wave strong-coupling regime we have observed other unusual facts already 



partially predicted in s-wave symmetry []lq-|22]| , but here enhanced by the specific properties 



of Bi 2212. Some of them can be discussed with reference to Fig. 2 where the real and 
imaginary parts of the renormalization function Z s (u) (Fig. 2a) and of the gap function 
A d (u) = V2A d (u) (Fig. 2b) are shown for T = 4.2, 99 and 140 K. In particular, one can 
observe that the function Real(Adiuj)) = Aai(u) goes to zero very abruptly for uj — > and 
T ~ T c (see the short dash curve of Fig. 2b), but, in general, at that temperature it is 
considerable different from zero at higher energies, whereas in the weak and intermediate 
coupling regime the whole curve goes to zero in uniform way at the increase of temperature. 
In addition, Fig. 1 shows that the normalized conductance is practically flat at T > T* but, 
at the same temperature, A^ (u>) is still different from zero, as it is suggested by dot and 
short dot curves of Fig. 2b that are calculated at T = 140 K. 

In the best fit case represented by the continuous line of Fig. 1 and at T = 4 
K, our theoretical model produces A = V2A (u — 0) = 21.22 meV that gives a ratio 
2A /k B T™ lc = 4.98, but the values of the gap edge A g = V2A g = v^real [A (u = A g (T))} 
and the peak Am are rather different: A g = 24.5 meV and Am = 25 meV. In Fig. 3a we 
can see the plot of A , A g and Am as function of temperature: the last one (solid triangles) 
does not go to zero at T c and, in addition, it grows considerably at larger T, being thus 
present also at T > T c . Particularly at T > T c /2 this peak is very different from the gap 
edge A g (solid squares). On the other hand, it can be seen in Fig. 3a that both A g and A 
have a BCS-like temperature dependency, but with a different T c of the order of 115 K and 



99 K, respectively. This effect is typical of a very strong-coupling regime while in the weak 
and intermediate coupling A g and A always coincide. 

Another very interesting and unusual phenomenon related to A g occurs for high values 



of A, |22|. When X s > 2, a T** exists for which at T > T** (in our case T** = 65 K) there 
is an energy value u*(T) that verifies the expression: real(Z s (u)*(T))) = imag(Z s (u)*(T))). 
In this case the quasiparticle approximation is no more valid and it is possible to have more 
than one solution of the equation A g (T) = real [A (uj = A g (T),T)}. In Fig. 3b we show the 
temperature dependencies of the two solutions A g i, A g2 of the previous equation, obtained 
for the spectral function of Fig. 1 and the parameters of the solid line curve of the same 
figure, together with the u*(T) curve and its approximate value u^ h (T) given by 

ul h (T) ~ [2vr/(l + A.)] / a 2 00 F(Q) • [f (fi) + n (fl)] dSl. (7) 

Jo 

This effect is also present in pure s-wave symmetry and is strongly enhanced by the particular 
shape of the spectral function of Bi 2212. In conventional strong-coupling superconductors 
the range of temperature where two solutions of the gap equation are present is of the order 
of 10~ 2 K [22], while here it is of the order of 50 K. Particularly in the present situation, 
this fact suggests that the gap edge is no more the most relevant physical quantity. As a 
consequence, the gap, experimentally determined from different measurements, is not the 
gap edge A g but an average over the frequency of the function A(u,T) with a weight 
factor depending on the property being measured ||22|| . This fact could have important 



consequences in the interpretation of the experimental data in Bi 2212. 

In summary, we have shown that the simplest approach to the pure <i-wave solution of 
the equations for the retarded strong electron-boson interaction reproduces very well and in 
the whole energy range the normalized conductance of optimally-doped Bi 2212, recently de- 
termined in our break-junction tunneling experiments 0,|J. The addition of a small amount 
of impurity scattering in the unitary limit accounts also for the small but finite value of the 
normalized conductance at zero bias |23| . In the best-fit cases, the calculated critical temper- 
atures differ no more than 6% from the experimental ones. The temperature dependency of 



the normalized conductance, calculated for the best-fit values of the electron-boson coupling 
constants, shows non-flat curves at T > T c up to T* f 40 K. This numerical result seems to 
suggest the presence of a strong-coupling "pseudogap" even in the DOS of optimally-doped 
Bi 2212 samples where, up to now, it has never been clearly observed. The evolution of the 
present work will be its extension to the non half-filling case, in order to compare the <i-wave 
solution of the Eliashberg equations to the DOS of underdoped and overdoped Bi 2212, and 
the comparison with other experimental data. Nevertheless, we believe the generality of 
the Eliashberg formalism and its independence of a specific assumption for the microscopic 
coupling mechanism gives the present results a hardly questionable validity. 
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FIGURE CAPTIONS 



Fig. 1 The theoretical normalized conductances for /^i = at various temperature values: 
T = 4, 60, 100, 120, 140 K and for fi* u = 0.0125 at T = 4 K, calculated with X s = 2 and 
Xd = 2.3. The optimally-doped Bi 2212 break-junction experimental data at T = 4.2 K are 
shown as open circles. In the inset, the a 2 F (^) Bi2 2i2 use d in the fit of experimental data is 
shown. 

Fig. 2 (a) The calculated values of ReZ s (uj) = Z s i(uj) and ImZ s {u) = Z s2 (uj) at three 
different temperatures T = 4.2, 99, 140 K; (b) the same as in (a) but for ReA d (uo) = A d i(uj) 
and ImA d (u) = A d2 (uj). 

Fig. 3 (a) Calculated values of A M (T) (peak of the DOS), A{u = 0,T) and A g (T) (gap 
edge) as function of temperature. The solid and dash lines represent the BCS dependence 
of the gap for T c = 115 K and T c = 99 K, respectively; (b) the two solutions of the gap 
equation A g i and A g2 , together with uj*{T) and its approximate analytical expression (short 
dash line) as function of temperature. 
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